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Introduction 


The effective population size (N,), usually computed 
from increase in inbreeding (AF), is a major criterion 
for classifying animal populations for degrees of 
(FAO 1998; Duchev et al. 
However, the interpretation of these parameters is 
their 
depends on completeness and depth of the pedigree 
of the studied population. And indirectly, for the 
effects of drift, mating policy, overlapping of genera- 
tions, selection, migration or population subdivision. 
Overall, these constraints and others move real ani- 
mal populations far from the theoretically idealized 
population referred to by genealogical tools (Fal- 


endangerment 


not straightforward because 


coner & Mackay 1996). 
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Summary 


We introduce a simple method to estimate effective population size from 
increase in coancestry (Ac) for all pairs of individuals j and k in a refer- 
ence subpopulation. An increase in pairwise coancestry for any pair of 
individuals j and k can be defined assuming that a hypothetical mating 
between them would give an individual with an inbreeding coefficient 
equal to cx, where cj is the coancestry coefficient between the individu- 
als j and k. The equivalent measure to discrete generations value (gjx) 
corresponding to the individual jk can be computed by averaging dis- 
crete equivalents generations of its parents (g; and gą). The mean 
increase in coancestry for all pairs of individuals in a reference subpopu- 
lation can be used to estimate a realized effective population size based 
on coancestries that would provide information on the effective size of a 
population under random mating. Performance of the new parameter 
was tested on simulated and empirical (horse) populations with different 
mating strategies and population structures. The routines needed to 
compute the introduced parameters have been included in a new ver- 
sion of the program ENDOG. 


To overcome these problems, Gutiérrez et al. 
(2008, 2009) have recently proposed the estimation 
of a realized effective size (N.; Cervantes et al. 2008) 
from the definition of the individual increase in 
inbreeding coefficients (Gonzalez-Recio et al. 2007). 

However, individual increase in inbreeding coeffi- 
cients are still affected by population structures 
(Gutiérrez et al. 2008). The computed effective popu- 
lation size using AF; accounts for all historical 
pedigree of the individuals and the obtained N, 
summarizes all the genealogical information of each 
individual. When structure exists (i.e. if local bottle- 
necks occur), AF; adjusts the recorded genealogies to 
account for this effect as soon as structure appears. 
This is a major issue in genealogical studies of real 
populations. 


2006). 


computation 
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It is well known in population genetics that under 
a regular system of inbreeding, both the increase in 
inbreeding and the increase in coancestry reach an 
identical asymptotic value with the only exception 
of scenarios of permanent subdivision (Falconer & 
Mackay 1996). However, these theoretical assump- 
tions do not work correctly in real populations 
because structuring, owing to geographical distance, 
selection or local founder effects, is the norm in real 
pedigrees with sufficient size and depth, particularly 
when the populations are under selective processes. 
This population structure is permanent but incom- 
plete with intensity varying over time because of the 
success of new selected families, improvements in 
reproduction and transport facilities or changes in 
the structure of herds. This structuring is very often 
cryptic and in subpopulations (i.e. individuals shar- 
ing a major proportion of genes of a selected family) 
cannot be straightforwardly identified. 

When dealing with structured populations, Ne can 
be estimated from increase in coancestry instead of 
increase in inbreeding (Caballero 1994). In contrast 
with increase in inbreeding, increase in coancestry is 
negligibly affected by population structure (Gutiérrez 
et al. 2008), which basically would reflect the effect 
of drift because of the finite size of the population 
(Caballero 1994). On the other hand, estimation of 
N, from increase in inbreeding could be preferred 
when the focus is on inbreeding depression which 
depends on the inbreeding coefficient and not on 
the coancestry coefficient. 

Here, we extend the methodology previously 
shown for the individual increase in inbreeding 
(Gutiérrez et al. 2008, 2009) to define an increase in 
pairwise coancestry (Ac) between individuals j and 
k with the aim of making it possible to estimate Ne 
from increase in coancestry in real populations with 
generations overlapping. 


Material and methods 
Theoretical background 


The individual increase in inbreeding described 
by Gutiérrez et al. (2008) with the posterior modifi- 
cation to account for avoidance of self-fertiliza- 
tion (Gutiérrez etal. 2009) is defined as AF; = 
1— “W1 -— F; with g; the discrete equivalent genera- 
tion computed for the pedigree of each individual as 
the sum over all known ancestors of the term (%4)”", 
where n is the number of generations separating the 
individual from each known ancestor (Maignel et al. 
1996; Boichard et al. 1997) and F; the inbreeding 
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coefficient of an individual i. This method takes 
advantage of the parameter equivalent to discrete 
generations (g, Maignel etal. 1996; Boichard et al. 
1997) to solve the problem that, in fact, overlapping 
generations are the normal scenario in real popula- 
tions. The parameter equivalent to discrete genera- 
tions was previously used to estimate effective 
population size by regression of inbreeding coeffi- 
cients over generations instead of regressing them 
over year of birth (Gutiérrez et al. 2003). This over- 
comes the problem of differences in pedigree depth 
in individuals born on a given date. However, 
regression-based methods to compute N, have been 
shown to be seriously affected by changes in the 
breeding policy, sometimes leading to nonsense neg- 
ative estimates of effective population size (Cervan- 
tes et al. 2008; Gutiérrez et al. 2008). 

The AF; values of the individuals belonging to a 
reference population can be averaged to give AF 
and, from this, a mean effective population size N,, 
called realized effective size by Cervantes et al. 
(2008), can be straightforwardly computed as 
N.= = and a standard error of N, can be further 
computed from the standard deviation of the 
individual increase in inbreeding and the square 
root of the reference subpopulation size as 


Oy, = 2N oar Fe (Cervantes etal. 2008) assuming 


N. independent individuals. The method by Gutiérrez 
etal. (2008, 2009) has been shown to indirectly 
account for most effects leading to increases in 
inbreeding higher than those expected for the ideal- 
ized population because they affect the individual 
inbreeding coefficients (Cervantes et al. 2008). 

In a similar way, an increase in pairwise coances- 
try for any pair of individuals j and k, belonging to 
different sexes, can be defined assuming that a 
hypothetical mating between them would give an 
individual with an inbreeding coefficient equal to cj, 
where c; is the coancestry coefficient between the 
individuals j and k. All coancestries including those 
between animals with the same sex can also be 
taken into account for a different use. The pedigree 
depth (gj) corresponding to the jk offspring can be 
computed averaging the equivalent complete genera- 
tions of its parents (g; and gg). Therefore, increase in 
coancestry between any pair of individuals j and k 
can be computed as, Acj = 1 — EAT — Ge where cg 
is the inbreeding of a descendent from both, and g; 
and gx are the discrete equivalent generation for the 
parents. This also takes into account the absence of 
self-fertilization (Gutiérrez et al. 2009). By averaging 
the increase in coancestry for all pairs of individuals 
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in a reference subpopulation, we can estimate a real- 
ized effective population size based on coancestries 
as Ne =-= that would provide information on 
the effective size of a population under random 
mating. 

As described by Cervantes et al. (2008), given that 
an increase in coancestry exists for each pair of 
animals in the reference subpopulation, a standard 
error of the realized N,. can be further computed 
from the standard deviation of these increases 
in coancestry (c4), and the square root of the 
effective size of the effective number of pairwise 
coancestries in the reference subpopulation as 


re . : : : 
ON, = 2N ofic In this expression, it is 


assumed that not all the coancestries are indepen- 
dent but only those computed in a population 
formed by N,, individuals. The reliability of this 
standard error, and also that for N, is still not 
absolute given that it is assumed that increase in 
pairwise coancestries are assumed uncorrelated, but 
it can be useful as far as it provides an idea of the 
variability and a comparative reference. 

We demonstrate this methodology on two ad hoc 
simulated populations and on several linked horse 
studbooks which can be studied separately and 
together as a single subdivided population. 


Simulated data 


Three scenarios were created, in all of them the 
founder population was formed by the same number 
of individuals of both sexes. A total of 100 males 
and 100 females were born for each generation and 
acted as parents under random mating, with no dif- 
ferential viability or fertility. No restriction on the 
number and sex of offspring was established. The 
following assumptions were considered: 

Simulation (i) The population consisted of a total 
of 100 males and 100 females born in each genera- 
tion and acted as parents for the following discrete 
generations. After 10 generations, the population 
was split in 4 different subpopulations consisting of 
25 males and 25 females evolving separately for 40 
additional generations. 

Simulation (ii) In this case, the initial population is 
subdivided from the foundation in 4 different sub- 
populations consisting of 25 males and 25 females 
each born in each generation and acted as parents 
for the following discrete generations, evolving sepa- 
rately until generation 10. After that, the four sub- 
populations were merged into a single population 
formed by 100 males and 100 females evolving 
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together under random mating over 40 additional 
generations. 

Simulation (iii) The population consisted of a total 
of 100 males and 100 females born at each time per- 
iod limiting the replacement of reproducing individ- 
uals to 50 males and 50 females at each time period 
and allowing the reproducing individuals to have 
offspring at two consecutive periods (overlapping 
generations). After 20 periods, the population was 
split into two subpopulations consisting of 50 males 
and 50 females evolving separately for 80 additional 
periods. 


Real data 


Six breeds previously analyzed by Cervantes et al. 
(2008, 2009) were employed to demonstrate the 
methodology using data from real populations, 
including parental and derived populations. The data 
included information from the studbooks of three 
pure horse breeds [Spanish Arab (SA), Spanish Pure- 
bred (SPB) and Spanish Thoroughbred (TB)] and 
three composite breeds [Spanish Anglo-Arab (dAA), 
Hispano-Arab (dHA) and Spanish Sport horse breeds 
(dSSH)] and collected between the foundation of the 
studbook and 2004: 

(i) The pedigree of the Spanish Purebred horse 
(SPB, Andalusian horse), totalling 140 629 individu- 
als (Cervantes et al. 2009). We also included the Car- 
thusian strain which is a subset of the pedigree of 
SPB with their own management structure, includ- 
ing the founder individuals classified as Carthusian 
because of their particular historical origin and the 
individuals resulting from matings between Carthu- 
sian parents (Valera et al. 2005). 

(ii) The Spanish Thoroughbred (TB) pedigree, with 
a total of 33 463 individuals (Cervantes et al. 2009). 
(iii) The Spanish Arab horse (SA) pedigree included 
18 880 individuals (Cervantes et al. 2009). 

(iv) The Anglo-Arab (dAA) breed that was formed 
by continuous crossing between Arab horses and 
Thoroughbred (TB) and totalled 8289 dAA animals 
(Cervantes et al. 2009). 

(v) The Hispano-Arab (dHA) breed that was formed 
by continuous crossing between Spanish Purebred 
(SPB; Andalusian) horses. The total number of dHA 
animals included in the dataset was 3394 (Cervantes 
et al. 2009). 

(vi) The Spanish Sport horse (dSSH), is a composite 
population that includes any individual participating 
in sports competitions that is not registered in any 
other studbook and, therefore, an individual belong- 
ing to SA, SPB, TB, dAA or dHA or any other breed 
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that can act as parent in the dSSH studbook (Cer- 
vantes etal. 2009). The total number of dSSH 
animals included in the dataset was 7099 (Cervantes 
et al. 2009). 

Because of the significant genealogical connected- 
ness among pedigrees analyzed by founder contribu- 
tion and by coancestry matrix (for more details see 
Cervantes et al. 2009), the six horse studbooks were 
further merged and also analyzed as a single set of 
populations. 

The horse pedigrees used here are expected to be 
useful to assess the ability of the introduced method- 
ologies to identify cryptic structures within popula- 
tions. From the history of the breeds, it is known 
that (i) within the SPB breed, the Carthusian popu- 
lation has a separate management, and within Car- 
thusians, some studs such as that of ‘Hierro del 
Bocado’ are particularly important in size and as a 
source of genes (Valera et al. 2005); (ii) the SA breed 
has a deep pedigree which is still open for Arab 
horses from foreign countries and, within this breed, 
the Military Stud of Jerez de la Frontera has a partic- 
ular breeding significance given that many farmers 
only select individuals from this origin to mate their 
animals (Maxwell 1995); (iii) the TB population has 
the shallowest pedigree of the parental horse popula- 
tions studied (Cervantes et al. 2009); (iv) the level of 
integration of the parental breeds into the derived 
breed (dAA, dHA and dSSH) is different and, as a 
whole, the set of horse populations analyzed here 
has unbalanced contributions from the different 
individual horse breeds. 


Reference subpopulations 


Changes in Ne, which are a consequence of varia- 
tions in the population structure or mating policies, 
can be addressed by defining several reference sub- 
populations across time periods and plotting the 
obtained values by time. However, if only the ani- 
mals born in the time period equivalent to the last 
equivalent generation interval are considered, the 
computed N, is likely to correspond to the present 
reproductive subpopulation. Thus, in simulations i) 
and ii) we used each generation as reference sub- 
populations (RS) to analyze the effective size trends, 
while in the simulation iii) individuals of the last 
two generations were used to obtain values of a 
reproductive subpopulation for subdivision and over- 
lapped generations combined scenario. In scenario i), 
the theoretical N, per discrete generation in terms of 
coancestry, excluding self-fertilization, is the number 
of individuals + % (200.5), and in terms of inbreed- 
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ing the same value before subdivision and 50 after 
that. Regarding scenario ii), the theoretical N, per 
discrete generation at the time of subdivision is 
200.5 in terms of coancestry and 50 in terms of 
inbreeding; after subdivision, theoretical N, both in 
terms of inbreeding and in terms of coancestry have 
the same value (200.5). Finally, for scenario iii) the 
theoretical N, in the last two generations in terms of 
coancestry can be derived from the expression 
Ne = Wet (Falconer & Mackay 1996), where Nc 
is the number of reproductive individuals included 
in the reference subpopulation (100), Vgm and Vey 
are, respectively, the variances of family sizes of 
reproducing males and females (Vy, = Vx, = 2 under 
random conditions), and L is the generation length 
in units of the specified time interval (1.5), Ne 
equals to 150. And in terms of inbreeding can be 
derived from the expression for unequal numbers of 
animals breeding in successive generations 


1 


A N TON (Falconer & Mackay 1996), where t is the 
number of considered generations and N; either 100 
or 50 depending on whether before or after subdivi- 
sion. In this case N, equals to 83.3. 

Regarding the real data, reference subpopulations 
were defined in real populations as the individuals 
born in the last 10 years from the available records 
for all the pedigrees. Pairwise coancestries between 
individuals were computed between all males and 
females in the RS except for SPB and the set of 
horse populations from which we selected a random 
sample of individuals because of computation time 
limitations. 

All the analyses were performed using the ENDOG 
program (current version v4.8) (Gutiérrez & Goy- 
ache 2005). 


Results 


Figure 1 shows the trend of realized effective popu- 
lation size estimated from AF; and from Acj in the 
simulated datasets (i and ii). In Figure la, both the 
effective sizes, as expected, are very similar until 
generation 10 where the simulated population was 
split into 4 subpopulations. After that the ratio Ne 
reached values over 3, and its trend showed that the 
real number of subpopulations can be reached 
asymptotically. In Figure 1b, the opposite scenario is 
shown: the ratio Ne value was nearly four until gen- 
eration 10 where population structure disappeared 
and both estimations lead to the same values. Com- 
mon founders who had old relationships seem to 


have some weight on later generations (la). Plot 1b 
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Figure 1 Trends of realized effective population size (N., squares 
and effective size based on an increase in pairwise coancestries 
(Nec circles) in the simulated example (i) (a; population split in 4 sub- 
population after 10 generations) and in the simulated example (ii) (b; 
population subdivided in four subpopulations until generation 10). The 
x-axis shows generation number and the y-axis the effective sizes 
values. 


shows that admixing unrelated populations immedi- 
ately gives a ratio similar to that corresponding to a 
single population. Regarding simulation iii) and 
using the last two periods as reference population, 
the estimated N, was 84.6+ 1.43 and N, was 
149.5 + 1.14, roughly the same values as the theo- 
retical effective sizes (83.3 and 150, respectively). In 
this case, the ratio between the parameters was 1.8. 

Table 1 gives the realized effective population sizes 
based on inbreeding and on coancestries for the RS 
defined in the real datasets analyzed. The values 
estimated for N, varied from 19.8 for the Carthusian 
horse strain to 332.1 for Spanish Thoroughbred, 
whereas the estimates for Ne were always higher, 
varying from 24.2 for the Carthusians to 819.1 for 
SSH. As expected, parental horse breeds reached 
lower estimates than composite breeds except for TB 
that had the shallower pedigree among the parental 
horse populations thus probably giving biased esti- 
mates of effective size. The estimates of N, had 
higher standard errors than Ne particularly for the 
composite horse populations, but remember that 
these are only standard errors. 

The ratio Ù% in SA and SPB was similar and 


Ne 
roughly 1.5, with the Carthusian horse value of 1.2 
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which is lower than that for SPB. In the Spanish 
Thoroughbred, the ratio was close to 1. Regarding 
the composite breeds, the lower ratio for dHA was 
1.1. This is considered as an unstructured popula- 
tion. However, a noticeable ratio of 4.6 was assessed 
for the dSSH population. Finally, the set of horse 
populations, including all the available breeds, had a 
ratio of 1.9. 


Discussion 


It is well known that many factors lead real popula- 
tions to increases in inbreeding higher than those 
expected in the idealized population as a conse- 
quence of preferential matings (Falconer & Mackay 
1996). However, increases in coancestry basically 
reflect the drift caused by the finite size of the popu- 
lation (Chesser et al. 1993; Caballero 1994; Wang & 
Caballero 1999). Thus, discrepancies between increases 
in inbreeding and in coancestry, except for mini- 
mum coancestry strategies, can be interpreted as 
cryptic population subdivision. 

Population structuring is one of the major factors 
leading to increases in inbreeding higher than that 
expected under random mating regarding population 
size (Fernandez ef al. 2008). On the other hand, 
minimum coancestry mating systems lead to lower 
increases in inbreeding than in coancestry that 
would be expected under random mating. In most 
real populations, mating policies, breeding goals or, 
simply, geographical distance prevent optimal mat- 
ings thus leading to population substructures. In 
such scenarios, non-random matings occur because 
reproduction is mainly carried out within subpopula- 
tions (Wang & Caballero 1999). Thus, substructuring 
greatly affects increase in inbreeding coefficients, 
while increase in coancestry remains approximately 
stable. When a permanent subdivision in a popula- 
tion is known, it is possible that this subdivision is 
taken into account and considered in the estimated 
parameters describing the population structure 
(Caballero 1994; Wang & Caballero 1999). 

However, in real scenarios it is not always easy to 
establish the presence of population structuring. 
Thus, for instance, geographical distances (and, 
indeed, other sorts of ‘distances’ such as type differ- 
ences) increase the probability that two individuals 
have common recent ancestors in the pedigree if the 
distance is short. Regarding herds, a population sub- 
structure leads to a situation in which some herds 
act as a nucleus, never receiving alleles from the rest 
of the population; these herds have higher increases 
in inbreeding than the others, while the mean 
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Table 1 The number of individuals born in the last 10 years of the available records for all the pedigrees (N:total, N: females, Nm: males), realized 
effective population size (N,) and the standard error (SE), effective size based on coancestries (N) and the standard error (SE), the mean of equiv- 
alent complete generations (t) and the standard error (SE), and number of computed coancestries (nc) for individuals included in the reference sub- 


population for each breed 


Population N Nm Ny Ne + SE Nec + SE t+ SE 

Spanish Purebred (SPB) 4460° 2209 2251 49.6 + 4.75 72.9 + 0.47 9.9 + 0.6 
Carthusian 1721 834 887 19.8 + 1.51 24.2 + 0.31 9.14 0.7 
Spanish Arab (SA) 6240 3032 3208 34.2 + 3.93 51.3 + 0.72 791.3 
Spanish Thoroughbred (TB) 3440 1910 1530 322.1 + 74.19 368.3 + 3.73 3.8 + 0.7 
Spanish Anglo-Arab (dAA) 2668 1297 1371 180.5 + 43.40 291.2 + 2.72 5.7 + 1.1 (2.5)? 
Hispano-Arab (dHA) 1783 925 858 128.5 + 42.66 139.6 + 1.50 7.7 + 2.0 (1.4)? 
Spanish Sport Horse (dSSH) 6055 2827 3228 176.7 + 71.68 819.1 + 8.93 3.4 + 2.0 (1.2)? 
Set of Horse populations 5472° 2679 2793 52.6 + 6.08 101.16 + 1.03 9.09 + 2.04 


Selected by random from a total of N = 89324. 


P Average equivalent complete generations computed excluding pedigree information of parents from parental studbook. 


“Selected by random from a total of N = 109510. 


coancestry grows equally among all of them (Vassal- 
lo et al. 1986). Elite groups within populations would 
also lead to the same situation; the mating between 
elite and non-elite animals would produce offspring 
that are excluded from the elite group. All these 
examples and many others provide higher increases 
in inbreeding than in coancestry as a consequence of 
population substructure. Under a regular system of 
inbreeding, and with a very deep knowledge of pedi- 
gree, both the increase in inbreeding and the 
increase in coancestry will attain an identical asymp- 
totic value except for the scenarios with permanent 
subdivision, but a long historical pedigree are not 
common in real livestock populations. 

If the aim is to analyze the inbreeding depression, 
the increase in inbreeding will provide a sound 
value, but when the goal is to ascertain the real 
effective number of breeding individuals in a struc- 
tured population, the estimation of the effective pop- 
ulation size from the increase in inbreeding performs 
poorly whatever the method used (Falconer & Mac- 
kay 1996; Gutiérrez et al. 2008). Inbreeding coeffi- 
cients are highly affected by population structuring 
and, in this case, the use of other methods such as 
that by Hill (1979), gathering the demography of a 
population, are more advisable (Cervantes et al. 
2008). Some other approaches developed to predict 
effective population size of structured populations 
are based on F-statistics (Wright 1969) and the vari- 
ances of offspring size of the identified subpopula- 
tions (see Wang & Caballero 1999 and Caballero & 
Toro 2002). These kinds of approaches, however, are 
not likely to be applied to real populations. 

In our study, we propose a simple method to 
estimate effective population size, directly, from the 
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increase in coancestry between each pair of individ- 
uals j and k (dcx) in a given reference population. 
This method is an extension of the methodology 
previously used for individual increase in inbreed- 
ing by Gutiérrez et al. (2008, 2009) that overcomes 
the difficulty of assigning each individual in a 
given population to a generation and, therefore, 
applicable to real pedigrees to obtain reliable esti- 
mates of effective population size based on coan- 
cestry (Ne). Moreover, a  between-individuals 
coancestry matrix, standardized for the complete- 
ness of pedigree, can be built and analyzed using 
other statistical tools, such as that based on princi- 
pal component analysis, to illustrate population 
structure (Jackson 1981). 

Furthermore, it has been shown that the compari- 
son between the new parameter and the individual 
increase in inbreeding gives information on the 
degree of population structure. N. and its brother 
parameter N, (Gutiérrez et al. 2008) are assumed to 
be measures of the same cumulated drift processes 
from the foundation of the population to the pres- 
ent. As they would be asymptotically equivalent in 
an idealized population, the disagreement between 
them mainly falls on their ability to assess the effect 
of preferential matings. In other words, the compari- 
son between N,. and N, would characterize the 
influence of preferential matings in the population. 
The method achieved reliable results in the three 
scenarios of simulated pedigrees giving the theoreti- 
cal value for each simulation. In the case of simula- 
tions i) and ii) (Figure 1), both Ne and N, reflect 
the situation produced by the subdivision effect. 
Moreover, when the subdivision is combined with 
a scenario overlapping generations, the estimated 
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Ne values approach the theoretical N, values better 
than those obtained using unadjusted methodol- 
ogies. 

The horse example is useful to discuss the infor- 
mation provided by both effective sizes. The avail- 
able Spanish horse populations can be studied 
separately or as a set of populations. Note that the 
Ne was computed in this case using coancestries 
between males and females to avoid computational 
limitations, but this procedure can also be useful to 
reflect a more realistic effective size from a predic- 
tive point of view if we are only interested in the 
next generation. The comparison between increases 
in inbreeding and in coancestry, allowed us to ascer- 
tain the presence of population structure within all 
the parental populations analyzed, except for TB. 
Consistently with the history of the populations, the 
ratio Ne/Ne significantly deviates from 1 thus 
reflecting the presence of a particular management 
within the Carthusian, SPB and SA populations 
because of the existence of diffused subpopulations 
originated by famous studs and preferred by some 
breeders (see Valera et al. 2005). This is more evi- 
dent in the case of SPB that includes in its studbook 
the Carthusian strain (Valera et al. 2005). The simi- 
lar values for both the effective sizes computed for 
TB should be interpreted with caution because they 
are more likely to reflect a shallow pedigree (only 
3.8 equivalent discrete generations) because of the 
continuous use of foreign TB stallions with unavail- 
able genealogies (Cervantes et al. 2009). 

The results from the composite horse breeds are 
more difficult to interpret. Both dAA and dHA had 
sufficiently deep pedigrees and assessing population 
structure from the last generation can give the idea 
of the unstructured breeds. The N,, was 1.6 higher 
than N, on the dAA pedigree, reflecting a similar 
structure to parental populations such as SPB or SA, 
but the standard error found was higher in this ped- 
igree. In dHA that had the deepest pedigree of the 
derived breeds, both values were very similar. This 
would characterize the composite dHA breed as a 
population with good pedigree information and ran- 
dom matings. However, the low effective size esti- 
mated by Cervantes et al. (2009) in the dHA breed 
suggests that this is more likely to be a consequence 
of the small proportion of matings available 
between individuals of the composite breed. In 
Table 1, we have included the equivalent discrete 
generations computed when parents from pure 
breeds are considered as founders to show that ped- 
igree within the composite breed is very shallow in 
this case. 
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Regarding the dSHH breed, we can assume that 
the recent foundation of its studbook and the lower 
pedigree depth (3.4 equivalent generations) would 
lead to increases in coancestry that are not translated 
into increases in inbreeding. However, within breed, 
individuals are selected for three different disciplines 
(Dressage, Show jumping and Eventing). These mul- 
tiple selection policies probably lead to different pref- 
erential matings with specific stallions to improve 
the desired performance given an overall increase in 
inbreeding with no parallel increase in coancestry 
(with Ne almost 5 times higher than N,). Therefore, 
original founder contributions seem to be kept by 
mating individuals always within a selection objec- 
tive thus maintaining a high equivalent structuring 
value. 

The results from the example of the set of horse 
populations gave a ratio between the effective sizes 
close to two. SPB individuals are nearly 90% of the 
total set. This extremely high SPB representation in 
the total set leads to a situation in which the other 
breeds appear as isolated individuals in a big and, 
more or less, connected population. 


Conclusions 


The estimation of effective size from increases in 
pairwise coancestry introduced here has been shown 
to perform reasonably well in real populations. This 
can be compared with other methods to compute N, 
from the increases in inbreeding, namely that from 
the individual increase in inbreeding to characterize 
the structuring degree. Population managers would 
prefer one or the other methods if descriptive or pre- 
dictive purposes are involved, but both of them pro- 
vide interesting information. Note that assuming 
some inbreeding, structuring would allow lower 
coancestry and consequently less drift. The compari- 
son between the increases in coancestry and in 
inbreeding has performed reasonably well in both 
the simulated and the real examples used in this 
study. These enable us to consider the population as 
a whole to identify whether the mating policy car- 
ried out in the population tends to create substruc- 
ture (lines) within it or not. 
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